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ABSTRACT 


The transmission of acoustic energy from a tapered 
fiuid layer inte an undertiyime layer of fluid was 
investigated. Pulsed transmissions at frequencies of 
75.85 kHz and 107.5 KHz were used to insonify the upper 
layer while a small probe transducer was used to measure 
Signal amplitude variation with, depth in the lower layer. 
The experimental data collected indicate that a directional 
transmission of acoustical energy into the lower layer 
occurs over regions corresponding to the individual cutoff 
depths for normal modes in the upper layer. A simple method 
Ger Gsvtimetine the position of these regions is discussed 


and the effects of frequency and upper layer slope variations 


ame reported. 
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I. ZUNPRODUCTION 


The propagation of sound in shallow water has been the 
subject of considerable theoretical and experimental research 
beginning with the study of propagation in layered liquids 
by C. L. Pekeris [Ref. 1]. In the Pekeris problem a point 
source and point receiver are positioned in a non- 
absorbing fluid layer of constant depth which overlies 
a semi-infinite layer of fluid having higher sound speed 
and density. The dominant portion of the solution posed 
by Pekeris is a summation of normal mode contributions with 
cylindrical spreading losses. Numerous studies of acoustic 
propagation in two-layered fluid models have subsequently 
been conducted to verify or improve upon the Pekeris solution. 
One example of such research is the experimental work 
reported by Knudsen [Ref. 2] which demonstrated that the 
theory generally predicated the observed behavior within 
experimental error. 

When the depth of the top layer of fluid is not constant, 
the Pekeris solution is no longer valid and the problem 
becomes considerably more complex. Williams and Lewis [Ref. 3] 
posed a solution which postulates that the normal modes adjust 
themselves adiabatically to the local depth when the depth 
varies slowly. This hypothesis was subsequently supported 
by the results of experimental work reported by Eby, Williams, 


Ryan, and Tamarkin [Ref. 4]. 





As described in Ref. 5 and elsewhere, each normal mode 
may be thought of as a pair of plane waves propagating in 
a duct or waveguide by successive reflections at the 
boundaries. These plane waves are incident upon the 
boundaries at an angle which is determined by the charac- 
teristic propagation vector of the corresponding mode. If 
the angle of incidence is greater than the critical angle 
for the given boundaries, propagation with relatively low 
attenuation will occur because virtually all of the acoustic 
energy is reflected. When the angle of incidence is equal 
to the critical angle, transmission of energy across the 
boundary begins and the corresponding mode is said to have 
reached cutoff. In the classic two-layered fluid problem 
with constant geometry, the directions of the mode propaga- 
tion vectors (and thus the angle of incidence) vary with 
the transmission frequency. It is therefore possible to 
define a cutoff frequency for the nth mode which is given 


by Knudsen [Ref. 1] and others as follows: 


recall) yp 


no Th 5 
‘i (c,/c5) 


jr. 5 


where h = depth of the upper layer 
ci ae= sound speed in the upper layer 
Cos sound speed in the lower layer. 


This means that at frequencies less than fa? the nth and 


higher numbered modes are evanescent. 





The foregoing equation can also be solved for a cutoff 
depth, ha? by rearranging variables and holding the 
frequency constant. This yields the minimum upper layer 
depth required to allow propagation in the nth and lower 
numbered modes. Based upon the observations reported in 
Ref. 4, this cutoff depth concept may be extended to deal 
with the problem of propagation in a slowly tapered upper 
fluid layer as suggested in Fig. 1. When the geometry 
of the tapered layer and the sound speed in both layers 
are known, it should be possible to predict, with reasonable 
accuracy, the regions in which significant transmission of 
acoustic energy into the lower layer occurs. 

The optical analog of the foregoing concept has been 
verified in an experiment reported by Tien and Martin [Ref. 6]. 
The waveguide in their experiment was formed by depositing a 
thin dielectric film with a tapered edge on a substrate 
having a higher refractive index. A light beam from a laser 
was coupled into the film and subsequent transmission into 
the substrate was observed. The transmission of light energy 
into the substrate commenced at a distance of several optical 
wavelengths in front of the cutoff region and was concentrated 
in a beam projected at an angle of 15 degrees below the film- 
substrate interface. 

The primary objectives of the experimental research 
reported herein are as follows: 

1. To verify that the transmission of significant 
acoustical energy into the lower layer of a two-fluid, tapered 


waveguide is confined to well-defined regions corresponding 





to the points at which the modal cutoff depths are reached 
in the upper layer. 

2. To determine whether or not the transmissions 
observed in these regions create a beam-shaped acoustical 
field in the lower layer. 

3. To determine some of the major effects produced by 
separately varying the transmission frequency and the slope 


of the upper layer taper. 
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II. APPARATUS AND EXPERIMENTAL PROCEDURE 


A. APPARATUS 

The experiment reported herein was performed in a 
specially constructed plywood tank at the Naval Postgraduate 
School. The tank was fabricated from 3/4-inch plywood with 
all interior surfaces heavily coated with Varathane (TM) 
itgutd plastic tO provide watertight integrity. The upper 
edge of the tank was rimmed with 2x4-inch stock and the 
sides were reinforced with 3/4-inch plywood struts to 
Provide Maer aicy. 

The determination of the shape and dimensions of the 
tank was subject to several considerations. Initial research 
indicated that the procurement of the desired fluids would be 
the major cost factor in this experiment. It was therefore 
decided to reduce the volume of fluids required by using a 
wedge-shaped tank rather than a conventional rectangular 
design. A second benefit of the wedge shape was that it 
helped to reduce interference due to acoustic reflection 
from the sidewalls of the tank. The length of the tank was 
chosen to accommodate the use of a tapered upper fluid layer 
Wetmeslepes varying from 0.02 to 0.05 over a minimum distance 
of one meter. The depth of the tank was chosen to provide 
for a lower fluid layer of at least 40 em so that it could 
be considered semi-infinite for pulsed transmission. The 


basic features of the tank are depicted in Fig. 2a. 
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In the typical shallow water problem the sloping bottom 
gives the water layer above it a tapered shape and the water's 
upper surface forms a level, pressure-release boundary. In 
a two-fluid model it is impossible to slope the bottom layer 
of fluid without separating the fluids with a rigid inter- 
mediate layer of glass or similar material. The top surface 
of the upper layer was therefore sloped by using a reinforced 
sheet of 3/4-inch plywood covered with foam-filled neoprene 
wetsuit material 1/8-inch thick. At the frequencies of 
interest, it closely approximates a pressure-release surface. 
A framework of aluminum angle beams was mounted on top of the 
tank and accurately leveled. The neoprene-covered board was 
then suspended from this framework and could be accurately 
positioned at the desired depths and angles. 

The acoustic source was constructed from a rectangular 
bar of barium titanate with the following dimensions; 
length = 2.7 cm., width = 0.9 cm., thickness = 1.28 cm.. 
Biased for resonance in its thickness mode, this source had 
a measured beamwidth of approximately 180 degrees in the 
vertical direction and /O degrees in the horizontal direction 
at a frequency of 75.85 kHz. A number of frequencies were 
Useole a the range from {5 "to 200 kHz. The signal’ to noise 
ratio at these frequencies was sufficient to avoid amplifica- 
Eron wand Tiltering difficulties. 

Tne receiver, a standard LC5-2 probe transducer manu- 


factured by CELESCO Industries, has a lead zirconate element 
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with a diameter of 0.093 inch mounted on the tip of a 
tapered stainless steel tube. The acoustic center of the 
element is located 0.125 inch from the tip of the probe. 

The transducer was mounted vertically in a 3/8-inch stainless 
steel tube bent at right angles to reach down the side of 
the tank and project outward toward its center. A pointer 
was positioned directly above the tip of the transducer 

by means of an aluminum rod attached to the steel tube as 
depicted in Fig. eb. The transducer mechanism was then 
attached to a carriage which slid on level rails mounted 
above the edge of the tank. This system allowed the trans- 
ducer to be positioned vertically with an accuracy of plus 
or minus one millimeter. The same accuracy was achieved in 
positioning the transducer horizontally along the centerline 
of the tank. All measurements were made using the aluminum 
framework previously described as a reference. 

The source transducer was driven by the amplified output 
of a Wavetek Model 116 voltage-controlled generator which 
was triggered by a General Radio Model 1310-B oscillator. 
The received signal was amplified bya factor of 40 dB and 
displayed on a multiple-beam oscilloscope. No filtering 
was necessary to maintain adequate signal to noise ratios. 

The lower fluid layer was a heavy brine with a depth 
of 43.5 em. The brine was produced by mixing 100 lbs of 
Sedium chloride (NaCl) and 150 lbs of calcium chloride 


(caCl,-2H,0) in approximately 120 gallons of tap water. The 
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3 and a sound 


resultant solution had a density of 1170 kg/m 
speed of 1740 m/sec at 23 °C. Olive oil was used to form 
an upper fluid layer with a maximum thickness of 4.0 cm 
(after the sloping board was positioned in the tank). 
Olive oil was chosen because it would not deteriorate the 


neoprene covering on the board and had relatively low density 


and sound speed: 905 ke /m> and 1425 m/sec, respectively. 


E. PROCEDURE 

The neoprene-covered board was positioned at an initial 
angle of 2.28 degrees relative to the surface of the brine 
layer and intersecting that surface at a distance of 28 cm 
from the wide end of the tank. (Refer to Fig. 3.) This 
formed a tapered duct in the upper layer with an initial 
depth of 4.0 em and a length of 100.4 cm. The source 
transducer was positioned at the narrow end of the tank 
with its vertical axis on the tank centerline and its 
horizontal axis at a depth of 2.0 cm in the upper layer. 
The distance from the face of the transducer to the 
beginning of the tapered section was 23.5 cm. 

The general experimental procedure was very straight- 
forward. The source was excited with a tone burst at a 
repetition rate of 100 cycles per second and pulse durations 
of approximately 100 microseconds. The probe transducer 
was then moved systematically to obtain the acoustic signal 
amplitude at selected depths and distances along the 


centerline of the tank. The pertinent components of the 
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received signal were identified by their variation with a 
combination of source and receiver movement when necessary. 
Horizontal distance measurements were referenced to the 

point at which the tapered section of the upper layer 
terminated (i.e., upper layer thickness equal to zero). 

The positive direction is taken when moving toward the source 
(upper layer thickness increasing). Vertical distance 
measurements were referenced to the fluid-fluid interface 


with the positive direction upward. 
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TET oo RE UES 


A. GENERAL OBSERVATIONS AND REMARKS 

The received signal was too complex to allow adequate 
resolution of its components when frequencies above 120 KHz 
were used. The vertical and horizontal beamwidths of the 
source transducer were considerably wider at the lower 
ferquencies ultimately utilized, but it was possible to 
identify the components of the received signal. Components 
of the signal related to paths involving single reflections 
from the sides of the tank, had arrival times that averaged 
0.06 msec later than the direct path arrival time. This 
was confirmed by computing multipath arrival times based on 
the geometry of the model. The relative amplitude of the 
side-reflected component was larger at positions corresponding 
to large source-receiver separations. The directionality 
of the source caused the reflected component to be relatively 
insignificant at shorter separation distances. The higher 
sound speed in the lower fluid layer allowed signal components 
related to a single bottom bounce path to arrive as much as 
0.6 msec earlier than the direct path component. The 
relative amplitude of this component was, however, 
considerably less than that of the side-reflected component 
due to the longer path it traversed and to transmission 
iosees at the fluid-fluid interface. Considerable care was 
taken throughout the experiment to reduce the effects of 


multipath interference by measuring the received pulse at a 
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time just prior to the arrival of the side-reflected 

component. Figures 4 and 5 depict samples of the received 
signal at frequencies of 75.85 kHz and 107.5 kHz, respectively. 
Some residual interference effects appeared to be present 

upon analysis of the data, but their effects were relatively 
small and did not seriously degrade the overall results. 

The source transducer was maintained at the same position 
during all data collection. Trial movements of the source 
for purposes of determining positional effects and received 
signal component identification were accomplished prior to 
actual data collection. The amplification factor for the 
received signal was also held constant at 40 dB so that all 
Signal amplitudes, V, refer to the uncorrected, peak-to-peak 
voltage measured on the oscilloscope. The estimated uncertainty 
of these measurements was +5 mv during the first stage of the 
experiment wherein signal amplitudes commonly exceeded 0.4 
volts. During subsequent stages, signal amplitudes were 
generally less than 0.2 volts and were measured with an 
e@ecuracy of £2.5 mv. The estimated uncertainty in the 
receiver position was plus or minus one millimeter in all 
directions throughout the experiment. Error symbols on 
Ene eragnical plots indicate the undertainty in the data. 

The lines drawn between data points were sketched as a 


visual aid only. 
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B. FE@RSTVSTAGE RESULTS 
The first stage of the experiment was conducted under 
the Gollowing conddtions: 


1. Upper layer depth before beginning taper .... 4.0 em 


2, Wistance from source Go beginning of taper ... 23.5 cm 
Se Siletiemet apices. ....u mw. wee oe... ew eee «=—60060398 

4, Length of tapered section ..... Re ines » 100.4 em 
ee er OWN WER. . . . GR. SS ee. me ee 36H. 85 CkKHZ 
CRE IE Lio SiO OMels ts Es 0s ees ee ees. ORL msec 
ae Relies repetition Tatew . . sews. . ste . Heer. soo ee . LOR) Az 


The received acoustic signal amplitude was measured at depth 
increments of 0.5 ecm along the longitudinal centerline of 
the tank, beginning at a position directly below the vertex 
of the upper layer taper (x = 0.0). These measurements were 
repeated at intervals of 5 cm along the centerline until a 
nerigenmtal distanceset 105 cmeizem the vertex wasereached. 
Additional measurements were taken at selected positions 
when better definitionwot the acoustic»fleld»was~-desired. 
All of the data collected during this stage of the experiment 
were obtained within a period of eight hours. 

The relative size of the receiver prevented a thorough 
sampling of signal amplitudes in the tapered section of 
the upper layer but the data obtained strongly indicate 
that first mode propagation predominated. The size of the 
source and its positioning at half-depth in the upper layer 
are factors that contribute to the suppression of the second 
and higher numbered modes so that the predominance of the 


first mode is reasonable. 
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A preliminary analysis of the data revealed that, at 
horizontal distances greater than 25 cm from the vertex, 
the decrease in signal amplitude was approximately exponential 
with depth. This was confirmed by making comparative plots 
of V (at Z = 0.0) multiplied by exp[-Z] versus depth (Z). 
Such behavior corresponds to the theoretical "tailing off" 
of the normal mode structure in the upper layer. 

A marked change in the variation of signal amplitude 
with depth in the lower layer occurred between horizontal 
distances of 25 cm and 15 cm from the vertex of the taper. 
Signal amplitude measurements at distances between 25 cm and 
20 em revealed no discernible pattern for the amplitude 
variation. Commencing at approximately x = 20 cm, however, 
the depth corresponding to maximum signal amplitude increased 
at a relatively uniform rate as the horizontal distance to 
the vertex of the taper decreased. The relative shapes of 
the plots of signal amplitude vs. depth in this region are 
also indicative of a beam-like directional transmission of 
acoustical energy originating at a horizontal distance of 
approximately 18 to 20 cm from the vertex of the taper. 
(See Fig. 6.) The calculated cutoff depth for the first 
mode at a frequency of 75.85 kHz is 0.82 ecm. This thickness 
in the upper layer corresponds moe horizontal digtance of 
20.6 em from the vertex of the taper. This distance relates 
remarkably well to the estimated origin of energy transmission 


based upon the experimental data. 
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Graphs 1 through 4 are plots which were designed to 
facilitate an estimation of the apparent angle of energy 
transmission in the vertical plane. These plots of 
experimental data were constructed using the following 
sealing eriteria: 

1. Signal amplitudes (V) were assumed to have suffered 
cylindrical spreading losses over the distance (r) from the 
estimated origin of the apparent beam to the point of measure- 
ment. (Refer to Fig. 1.) The horizontal coordinate of a 
datum therefore represents the signal amplitude extrapolated 
back to a distance of one centimeter from the apparent origin. 

2. The depth (Z) at which the amplitude measurement was 
made was divided by the horizontal distance (x') from the 
assumed origin of the beam. The vertical coordinate of a 
datum therefore represents the tangent of a depression angle 
measured from the horizontal plane. 

3. The origin of the apparent beam was assumed to be 
at a horizontal distance of 20 em from the vertex of the 
taper. Transmission of acoustical energy into the lower 
layer may have occurred over a distance of approximately 
4.0 em centered about this origin. 

The assumption of cylindrical spreading was considered 
consistent with the spreading exhibited by the main beam 
of a line-like source at relatively short ranges. 

A sequential examination of Graphs 1 through 3 indicates 
that the beam-like transmission pattern was not well defined 


until a distance of approximately 10 cm from the estimated 
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origin was reached. At longer distances the pattern was 
better defined and the data plots were relatively similar 
when the possible variance due to measurement errors was 
considered. The finer structure of the plots is believed 
to be due, in part, to residual multi-path interference 
effects discussed earlier in this report. 

Because of the geometric factors involved in the scaling 
of the data, uncertainty in the estimated position of the 
beam origin can result in considerable error when deter- 
mining the values of x' and r. The relative magnitude of 
these errors increases as the horizontal distance from the 
estimated origin decreases. Data corresponding to values 
of x' less than 10 cm were therefore excluded when estimating 
apparent angles of transmission and beamwidths. These 
parameters were estimated by analyzing an overlay of the 
data corresponding to values of x' greater than or equal 
to 10 em, as depicted on Graph 4. A datum corresponding to 
the average maximum amplitude of the combined data was 
estimated by averaging numerically the individual maximum 


1/2) data. The vertical coordinate of the 


amplitude (Vr 
estimated datum was then used to provide the first estimate 
of the beam depression angle. The 3 dB down points were 
subsequently estimated in a similar manner using the 
horizontal coordinate of the maximum amplitude datum as a 


reference. The depression angles corresponding to these 


points thus provided an estimate of the apparent beamwidth. 
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A second estimate of the beam depression angle was obtained by 
bisecting the angle between the 3 dB down points. Table I 
provides a list of the estimated depression angles and 
beamwidths for all stages of the experiment. 

The foregoing estimates are, of course, subject to some 
error. Analysis of the geometry suggests that the apparent 
angle of transmission could be in error by about +2 degrees 
if the estimated origin of the transmission had an uncertainty 
of +2.0 cm. The estimated angular concentration of trans- 
mitted energy (beamwidth) can reasonably be expected to 
reflect a similar percentage of error. 

A period of five days elapsed between the conclusion of 
the first stage of the experiment and the commencement of 
further experimentation. During this period the temperature 
of the building in which the test tank was located was 
considerably lower than normal. A thin film of unidentified 
Substance formed along the oil-brine interface during this 
period, possibly as a result of the temperature reduction. 

A repetition of the first-stage amplitude measurements at 
selected positions revealed that the transmission patterns 
were relatively unchanged but signal amplitudes in the 
apparent beam were reduced by approximately 60%. Signal 
amplitudes in the upper layer at short distances from the 
source, however, were virtually unchanged. Graphs 5 and 6 
depict typical data obtained during a repetition of the first 


Stage. The data obtained at corresponding positions during 
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the first stage were reduced by 60% and superimposed on the 
graphs for comparative purposes using plus signs (+). 
Visual examination suggests that while amplitudes have teen 
reduced, the basic features of the transmission remain 
essentially umchanged. The resultant estimates of 
depression angle and beamwidth listed in Table I further 


support this contention. 


C. EFFECTS OF FREQUENCY AND SLOPE VARIATIONS 


@ae second stage of the experiment was carried. out 


” 
ae 


}4 


@ manner identical to the first stage with the following 


changes in conditicns: 


Re SC ec ew cw ttt ee 107.5 kHz 


ae. Pulse duration eeee#ee+eeeeeeee#e?eee#e ee @ OTs msec 
In accordance with the definition of the modal cutoff depth 


given previously, the following cutoff decths and distances 


were computed: 


h, = O«57S cm Sq = Gl.52 om 
A, = 1.734 cm . 3 42.55 em 
Ae = 2.89 em X2 = 72.6 cm 


Data were then collected at appropriate locations to 
determine whether beam-like transmissicns originated at 
Ghese positions. This resulted in the identification of 


transmissions related to the first and third modes. The 
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absence of the second mode was expected but the relatively 
large magnitude of the third mode transmission was somewhat 
supprising. 

Scaled plots of the experimental data were constructed 
in the same manner as those for the first stage. Graphs 7 
through 10 relate fhe to the first eae transmission. The 
origin of the apparent beam was estimated to be at a 
horizontal distance of 12 cm (+2 em) from the vertex of the 
taper which is fairly consistent with the theoretical 
distance. Refer to Table I for estimates of beam parameters. 
Graphs 11 through 13 relate to the third mode transmission. 
The estimated origin of the beam was at a horizontal 
distance of 70 ecm (+2 em) from the vertex mien is, again, 
consistent with the predicted position. Corresponding beam- 
width and depression angle estimates are listed in Table I. 

The final stage of the experiment was conducted under 


pace rollowing conditions: 


i Upper Wayerecdepthneprior CO Caper ....cssseeccccs 4.0 cm 

2 Distance from source to beginning of taper ..... 47.4 em 
SRE SMT AIMS! . 5 sn sc ove wes 0 6 Wlaleiete eee ce weceseee Oe052 

Wo Hemeth of tapered section .......ese00. so 0e ot see ORE Tn 
DEI WN I, occ cc hs od 6s bh See ce ee aes wesw scene Pow Oo ema 
Ge PUISe dUrat1OnN .....s08% Se a DAS ar ees Aas once ote SEONG »+. O.1 msec 
i Diese ereme Git OM Babe)... 6i6 4% 666 esse Se oa arene ete cee OOM Ea 


The procedures utilized were identical to those in the 


previous stages. 
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The cutoff depth for the first mode was the same as that 

computed for the first stage since the frequency was the: 

same. In this case, however, the distance from the vertex 

of the taper is reduced to 15.6 cm by virtue of the change 

in the slope. Based upon the data obtained, the estimated 
origin of the apparent beam was approximately 15 cm (+2 ecm) 
from the vertex. Graphs 14 through 16 depict the corresponding 
scaled data plots for the new slope conditions which were 
analyzed as before. The resultant estimates of beam width 


and depression angle are listed in Table I. 
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IV. CONCLUSIONS 


The data obtained during this preliminary investigation 
Were sufficient to support the following conclusions: 

1. The transmission of acoustical energy from a 
tapered fluid layer overlying an essentially semi-infinite 
fluid layer with higher sound speed is concentrated within 
relatively well defined regions along the fluid-fluid 
ceri ace . 

2. lhewmeegrvons of Concentrated transmission correspond 
to the modal cutoff depths which are frequency dependent. 
The positions of these regions may therefore be predicted 
with reasonable accuracy when the properties of the fluids 
and the slope of the taper are known. 

3. The transmission of acoustical energy into the lower 
fluid layer within fnese regions exhibits directional 
properties which indicate that a beam-forming process may 
occur over the length of the transmission region. 

4, Additional data for larger ranges of frequency and 
slope would be required to draw any firm conclusions 
regarding their influence on the beamwidth/depression angle 


of the beam sent into the bottom. 
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TABLE I 


Estimated 
F(kHz) Mode Slope Soe ie Teco, saa ae 
75.85 a 0.0398 ay 132 12355 1-4 
75.85% i 0.0398 18° 11.30 ease 5,6 
Of. a 0.0398 14.4° 1053" Ley 7-10 
OT <D 3 0.0398 Hoel ro" ser 11-13 
75.85 a O.052 16n 3° Ne SoS 14-16 


ee LL 


Method A: Using average maximum amplitude datum. 
Method B: Using-bisector of average -3 dB data. 


*Data from the repetition of the first stage 


(Estimated errors in all angles about +2°.) 
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